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The wandering exponent v for an isotropic polymer is predicted remarkably well by a 
simple argument due to Flory. By considering oriented polymers living in a one-parameter 
family of background tangent fields, we are able to relate the wandering exponent to the 
exponent in the background field through an e-expansion. We then choose the background 
field to have the same correlations as the individual polymer, thus self-consistently solving 
for v. We find v = S/(d + 2) for d < 4 and v = 1/2 for d > 4, which is exactly the Flory 
result. 



2 April 1993; Revised 2 June 1993 



email: kamien@guinness.ias.edu 



1. Introduction 



The theory of polymer conformations is at once elegant and confounding. Flory theory, 
which is based on dimensional analysis, does a remarkable job in predicting the wandering 
exponent v defined by ( [r(L) —r(0)] 2 } ~ L 2v [1]. Nonetheless, controlled approximations 
which attempt to improve on Flory theory for polymers are not nearly as good. Flory 
theory predicts v F = 1/2 for d > 4 and v F = 3/(d+ 2) for d < 4, which is exact in and 
above the critical dimension d c = 4, where the polymer acts as a random walk [2]. It 
has been shown [3] that v F is exact in d = 2 dimensions as well. In three-dimensions the 
resummed e-expansion gives v = 0.5880 ± 0.0010 [4], Flory theory gives v F = 0.6, and the 
most recent numerical simulations cannot distinguish between the two [5]. 

In this note we rederive the Flory exponent. We do this by considering oriented poly- 
mers interacting with a background directing field. Flory theory would predict the same 
value of v for oriented or non-oriented polymers. For a family of directing fields, parame- 
terized by A, we can derive an exact result for v in terms of A. Unlike other self-consistent 
analyses, we find the self-consistent exponent by matching exponents, not correlation func- 
tion prefactors [6]. We then choose A so that the tangent-tangent correlation function of 
the polymer scales with the same exponent as the background field correlation, thus self- 
consistently choosing A. We find the Flory value, namely v = 1/z = S/(d + 2), where z 
is the dynamical exponent, as we discuss below. It is amusing that the Flory result comes 
from an e expansion which happens to be, in this case, exact. This may suggest why Flory 
theory is so good. Additionally, it suggests how one might study tethered surfaces, where 
Flory theory is not so good and e-expansions are not so easy. 



2. Formulation 

We consider a d dimensional system with a background tangent field u, 



ZJr, s]=J\f 



/"> ,exp -^r ds '(^- Au(rM - s ' ) ) 



(2.1) 



r(0)=0 

where s labels the monomer along the polymer. Then the annealed partition function is 



^an(r, s) = J [du]Z u (r, s)P[u}. (2.2) 



A factor of A has been introduced to help to organize the perturbation expansion, and 
will, in the end, be set to unity. In addition, it is needed to make the units correct, i.e. 
[A] = LS- 1 . 

The vector field u(x,s) is a function of both space and the monomer label s, as if 
each monomer interacts with a different vector field, even though they could be near each 
other in space. Note that in mean field theory dr(s)/ds — Xu(r(s), s). Consider the 
tangent correlations of two monomers at two nearby points in space, A and B. If the 
polymer takes a short route from A to £>, then the tangent vectors along that length will 
be strongly correlated, since they must mostly lie along the line connecting the two points. 
On the other hand, if the polymer takes a long circuitous path while getting from A to 
£>, the two tangent vectors will not be very correlated. Since u(x, s) is the self-consistent 
tangent field, it must reflect this simple geometric argument. If u did not depend on s, 
the tangent-tangent correlations along the polymer would only depend on their distance 
in space, and would not respect the constraint relating the path to the tangent vectors. 

We now view the functional integral over r(s) as a quantum mechanical propagator 
in imaginary time. That is, Z u will satisfy a Fokker-Planck equation with the initial 
condition r(0) = 0. In this case the Euclidean Lagrangian is just our Hamiltonian. The 
Euclidean Hamiltonian is the Legendre transform of the Lagrangian. However, as with the 
theory of a point particle in an electromagnetic field, there is an ordering ambiguity. The 
functional integral produces the symmetric Hamiltonian. Because, in statistical mechanics, 
the Hamiltonian comes from a transfer matrix, we must take the symmetric ordering of the 
momentum and the velocity field, as is done in electrodynamics of a quantum mechanical 
particle [7]. 

The Fokker-Planck (Euclidean Schrodinger) equation obtained is: 

8Zu q s ' S) = DV 2 Z u (r, s) - h {V- [«(r, s)Z u (r, s)} + «(r, s) ■ VZ u (r, s)} . (2.3) 
If we write ip(r, s) = 9(s)Z u (r, s) and Fourier transform in space and time, we get: 

(-zu + Dk 2 ^u) = Mk) ~ l^X) J^J ^Uiik-qXk' + tfMM) (2.4) 

where u is the fourier variable conjugate to s, k is the fourier variable conjugate to x and 
ipo is the fourier transform of the initial conditions. The boundary condition r(0) = 
corresponds to ip = 1. Equation (2.4) can be solved recursively in powers of A. 
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Now we must consider the random field it. If u is to describe the self-consistent 
background, then, if the polymers have no ends (by either being cyclic or spanning the 
system), V-u = 0. With this constraint, we take 

( u t (k, u) Uj (k>, u') ) = + u') S d (k + k>) (2.5) 

We will self-consistently choose A at the end of the calculation. We have chosen the uj 2 
dependence in the denominator for simplicity. If we had chosen a dependence other than 
quadratic we could always capture the same relative scaling of u> and k by an appropriate 
choice of A in (2.5). 

We note an important simplification due to the constraint V • u = 0, [8] namely 
that there is no difference between taking u to be quenched or annealed. The quenched 
probability distribution for r is just: 

^■"'/w /S,/ 1 ' 1 (2 ' 6) 

However, by integrating (2.3) over space, we find that 

^ J d d rZ u (r,s) = J d d rV-[DVZ u (r,s)-Xu(r)Z u (r,s)]. (2.7) 

where we have used the fact that V-it = 0. By gauss's law, the integral on the right hand 
side will vanish since Z u (r,s) will fall to at infinity. Thus the normalization of Z u is 
s-independent. Since J d d rZ u {r 1 0) is a constant, independent of u, so is J d d rZ u {r, s), and 
it factors out of the functional integrand in (2.6). But then we see that 

P (r J[du}Z u (r,s)P[u] 
-ran (r, s ) 



f d d r J[du]Z u (r, s)P[u] 

J[du]Z u (r, s)P[u] (2.8) 
Jd d rZ u (r, s)J[du]P[u] 

= P qu (r, s). 

Because of the directed nature of the propagator, it is possible to argue that z* = 2 — e 
to all orders. Since this problem can be viewed as quenched disorder, u will not suffer any 
nontrivial rescalings, and so A and B will only rescale trivially. In addition, as we will 
argue in Appendix A, A does not get renormalized at any order of perturbation theory. 



3. Perturbation Theory and the Renormalization Group 

We analyze this model within the context of the dynamical renormalization group, 
along the lines of [9]. For a renormalization group with parameter £, we rescale lengths 
according to k' = e t k and uj' = e 7 ^u;, where 7(f) is an arbitrary function of £ to be 
determined later. When integrating out high-momentum modes in a momentum shell 
Ae~ £ < k < A, we integrate over all values of uj. We can choose the field u to have 
dimension 0. We find differential recursion relations: 

^) =BW (-2 + Z+ ^(l-> <3,«) 

d\{t) 



di 



= X(£)(-l + z) (3.16) 
^P-=B(£)(d-2A + z) (3.1c) 

where z(£) = Y(£) and Ad = 2(4-7r) d / 2 /T(d/2) is a geometrical factor. Putting these 
together, we can get a recursion relation for g = A 2 (l — ^)A d /2BD : 

| =»(«-») (3.2) 

where e = 2A—d. Thus in d > 2A dimensions there is a stable fixed point at g = 0, whereas 
if d < 2A then there is a stable fixed point at g = e. In the details of the calculation, 
it is essential that A < 2. Succinctly put, this is because the pole in uj which is used to 



evaluate the self-energy correction (a; = a/ (B/A)k A ) must dominate the diffusive part of 
the propagator (Dk 2 ) at small k. We will check that this constraint holds when finding 
the self-consistent value of A. 



4. Self-Consistent Value of A 

Choosing D to be fixed simplifies calculations (though does not change the results), 
so we choose z{£) = 2 — g{£). As we show in Appendix A, at the nontrivial fixed point, 
z = 2 — e exactly. Using this we find the following two scaling relations for the position 
and velocity correlations: 

([r(L)-r(0)] 2 > = 2dDL 2 / z (4.1) 



and 

dr(L) dr(0) AdD(2-z) (2/z) _ 2 



Now we would like to look at the tangent field correlations, corresponding to (4.1) 
and (4.2). We would like to evaluate 

f d d k doj e ik ' r ^e~ iujs (d — 1) 
(Mi( r W , 3 K(0,0)> = y^ F - Aw2+B [ 2A («) 

Since (r(s)) = 0, we take r(s) = in (4.3) and then check that the corrections are 
themselves consistent. We find 



= (d-l)A d Y(l-e)A^I^ L{A _ d)/A 



< uMs),s) Ul (0, 0) ) = ^ ' a wd/2AA ^ (A " dj/A (4-4) 



Matching the scaling exponents, we have 

A - d 2-2z 4A - 2d - 2 



2 - 2A + d 



(4.5) 



Solving for A, we find A = d/2 or A = (d+2)/3. The former gives e = 2(d/2) -d = 0, and 
hence z = 2. In this case (4.2) is incorrect, since the second derivative of (4.1) vanishes. 
Thus, we must choose A = (d + 2)/3. We note that this solution satisfies both A < 2 
and 2A > d for d < 4, and so we can say that the critical dimension of this model is 
d c = 4. Above this, e is negative, and we return to the simple random walk fixed point, 
z = 2. Finally, we compute z = 2 — 2 A + d = (d+ 2)/3, in complete agreement with Flory 
theory. We point out that the matching of exponents breaks down when d = 4, for if (4.1) 
is described solely by an exponent, without logarithms, there is no matching to do - that 
is, the exponent in (4.2) would be 0. Thus, we do not expect, and in fact do not, reproduce 
the correct logarithmic corrections to scaling. 

If we expand the complex exponent in the integrand of (4.3) in powers of r(s) we 
would find an expansion of the form 

gMK^M («) 

where we have suppressed the indices. This sum only contains even powers since the 
integration over k will eliminate odd powers of k by rotational invariance. Since ( r(s) 2 ) J ~ 
s 2j / 2 and each power of k 2 will produce a factor of s _2 / A , we find that the higher order 
corrections scale the same way as the leading term if A = z, which, in fact, it does. Thus 
the approximation is truly self-consistent. 



5. Conclusions 



One might imagine adding to this model in a number of ways. One possibility is to 
add an explicit self-avoiding term to (2.1). This would only cause the same complications 
present in Flory theory. Another possibility would be to add a small divergence to the 
field u, corresponding to polymer heads and tails [10]. We would then have 

( u) Uj (k', J) ) = ^"f^f + + k') (5-1) 

where a is close to, but not equal to 1. In this case we find that new graphs arise which 
spoil the exact argument in Appendix A. Moreover, some of these new graphs will diverge 
logarithmically when A = Ap = (d + 2)/3. One might add a small correction 5 to Ap and 
do a 5 expansion around 6 = 0. Unfortunately, now, there is a new fixed point, and the 
self-consistent correction 5 is not small. 

It would have been more natural to consider a walk in a random potential and then 
self-consistently choose the two-point correlation of the potential to match the density- 
density correlation. It is, unfortunately, notoriously difficult to study polymers in a random 
potential in an arbitrary dimension [11]. Here we have exploited the solubility of a random- 
walker in a random velocity field. 

It is perhaps a curiosity that self-avoidance was not involved in this calculation. In 
fact, by mapping the system to quantum mechanics, we have actually mapped the system 
to a directed walk in an external vector field, ignoring the energy cost of self-intersections 
entirely. The self-consistent vector field is a strange object, as it depends not only on the 
polymer position r(s), but also the point along the polymer s. This is necessary from a 
geometric point of view, and suppresses tangent vector correlations between two monomers 
far apart along the polymer sequence. In (2.2) we can imagine integrating out the velocity 
field u. The details of this are in Appendix B, and the resulting free energy looks similar 
to that for a self-avoiding walk, with an energy cost for self-intersections. Reproducing the 
Flory exponent then, suggests that the Flory theory may be more robust, and that the 
model studied here is in the same universality class. The quality of the Flory prediction 
may lie in the fact that the self-consistent analysis here incorporated an exact result, within 
the epsilon expansion. 
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Appendix A. Exact Value of z 

In the above discussion we used the result that z = 2 — e at the fixed point, to all 
orders in e [12]. This result is similar to that in [9] where Galilean invariance assured that 
the scaling field of the interaction would only rescale by its naive dimension. Our argument 
will be perturbative. The first requirement is that u only rescales by its naive dimension. 
Since the quenched and annealed problems are the same, we can regard uasa quenched 
random field. While it is typical to absorb the coupling A into u, we will not do so here, 
and hence it will not be non-trivially renormalized. We also note that (2.3) is linear in ip 
so any nontrivial rescaling of ip will not affect A. 

We consider a general graph with an incoming ip line, carrying momentum p, and 
outgoing ifj* line, carrying momentum p' , and an outgoing u line, carrying momentum 
p — p'. The incoming ijj line must first meet a vertex with u. At this vertex, let u carry 
away momentum k\. The contribution to the graph from this vertex is then proportional 
to: 

P„ {(h) 2 5^ - fcffc?) = Vil {(h) 2 5^ - k»k{) (A.l) 

Thus the graph will be explicitly proportional to p M . Similarly, the outgoing if;* line will 
emerge from such a vertex. If the u line carries momentum k 2 , then it will contribute a 
term proportional to: 

(p' - k 2 ) v {{k 2 ) 2 5™ - k^kl) = (p'Y ((M V CT - k^kl) (A.2) 

due to the transverse nature of the u propagator. Thus the graph will generate a term 
with at least two powers of the external momentum. This will not then renormalize A since 
it is the coefficient of a term with only one power of momentum. Indeed, this graph will 
generate terms, which by power counting, are irrelevant operators. Hence, the recursion 
relation for A will simply be 
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to all orders. 
Now, 



^l = B(d-2A + z) (A.4) 
and if D has the recursion relation: 

^=D(-2 + z + f(g)) (A.5) 

where f(g) represents the perturbation expansion renormalizing D, then g must have the 
recursion relation: 

^ = g [2{-l + z}-{d-2A + z}-{-2 + z + f(g)}} 

= g[2A-d-f(g)} ( A -6) 

= g[e-f(g)} 

since the dimension of g ~ X 2 /BD will include a contribution from S and .D. Thus, if 
there is a fixed point, then f(g) = e to all orders in g. Hence, we see that at the fixed 
point z = 2 — e to all orders. 



Appendix B. Effective Self- Avoidance Interaction 

Starting with (2.2), the coupling of the polymer to the random field can be written 

*int = ^ / ***** ^(*) - s ) • ^ (B - 1} 

Upon integrating out u, we find to order A 2 (note that this will not include the 0(X 2 ) term 
in (2.1)) 

l2 f d d k , , . 



int " / (2yr) d 



-v^* A l— 'I , 2 fc^^-rrfW^ dri^dr,^) (B ' 2) 



16Dv / ABfc A + 2 V J J/ ds ds> 

Because of phase oscillations, the integral over k is small if r(s) — r(s') is large. If r(s) — 
r(s') = we can replace k 2 5ij by kikjd in the k integration, through rotational invariance. 



Making this substitution then is a good approximation, and the corrections are suppressed 
due to the oscillations. We then have 



int 



16DVAB J (2tt) 
\ 2 {d-l) f d d k 



X 2 (d-1) f d d k 




laD^AB J (27r) d ' 
X 2 (d-1)VB f d d k 




(B.3) 



16DVA* J (2tt) 



dsds' k A ~ 2 e ik - [r{s) ~ r{ - s ' )] e~^^ )k ^\ s ~ s '\ 



+ constant 



Within the self-consistent approximation we found A = (d + 2)/3, so A — 2 = d — 2 A = — e. 
We now expand in powers of \s — s'\, since by self-avoidance large values of \s — s'\ will 
typically be accompanied by large values of r(s) —r(s') which will suppress the integral. In 
addition we can expand in powers of e assuming that away from the critical dimension d = 4 
we only have a slightly modified potential, corresponding to a renormalized interaction. 
The first term in a double expansion in powers of e and \s — s'\ is 



looking very similar indeed to a self-avoidance term. We can also see that the repulsion is 
proportional to A 2 , which itself is proportional to the expansion parameter of the model 
discussed in this paper. 




(B.4) 
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